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Abstract 

The (axis-parallel) stabbing number of a given set of line segments is the maximum number 
of segments that can be intersected by any one (axis-parallel) line. This paper deals with finding 
perfect matchings, spanning trees, or triangulations of minimum stabbing number for a given 
set of vertices. The complexity of finding a spanning tree of minimum stabbing number is one 
of the original 30 questions on "The Open Problems Project" list of outstanding problems in 
computational geometry by Demaine, Mitchell, and O'Rourke. 

We show A/^T'-hardness of stabbing problems by means of a general proof technique. For 
matchings, this also implies a non-trivial lower bound on the approximability. On the positive side, 
we propose a cut based integer programming formulation for minimizing the stabbing number of 
matchings and spanning trees. From the corresponding linear programming relaxation we obtain 
polynomial-time lower bounds and show that there always is an optimal fractional solution that 
contains an edge of at least constant weight. We conjecture that the resulting iterated rounding 
scheme constitutes a constant-factor approximation algorithm. 

ACM Classification: F.2.2 Nonnumerical Algorithms and Problems. 
AMS Classification: 68Q17, 68U05, 90C27. 

Keywords: Stabbing number, matching, spanning tree, triangulation, complexity, linear program- 
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1 Introduction 



Objective Functions. Many problems in combinatorial optimization, algorithmic graph theory, or 
computational geometry deal with minimizing the length of a desired structure: given a set of vertices, 
find a set of line segments of small total length, such that a certain structural condition is maintained. 
Among the most popular structures are spanning trees, perfect matchings, or (in a planar geometric 
setting) triangulations of minimum total length. Other geometric problems give rise to other objective 
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functions: for example, one can ask for the total turn cost between adjacent line segments; e.g., see 
0. 

When dealing with structural or algorithmic properties, another possible objective function is the 
stabbing number, for a given set of line segments, this is the maximum number of segments that are 
encountered (in their interior or at an endpoint) by any line. If we consider only axis-parallel lines, 
we get the axis-parallel stabbing number. A closely related measure defined by Matousek |[20l is 
the crossing number, which is the number of connected components of the intersection of a line with 
the union of line segment^ If there are no connected components of collinear segments (which is 
the case for matchings), the crossing number coincides with the stabbing number. When considering 
structures like triangulations, the crossing number is precisely one more than the maximum number 
of triangles intersected by any one line. 

Stabbing problems have been considered for several years. The complexity of many algorithms 
in computational geometry is directly dependent on the complexity of ray shooting; as described by 
Agarwal H], the latter can be improved by making use of spanning trees of low stabbing number. We 
will sketch some related results further down. Most previous work on stabbing and crossing problems 
has focused on extremal properties, and little has been known about the computational complexity 
of actually finding structures of low stabbing number, or possible approximation algorithms. In fact, 
settling the complexity of Minimum Stabbing Number for spanning trees has been one of the original 
30 outstanding open problems of computational geometry on the list by Mitchell and O'Rourke ||2TI . 
(An up-to-date list is maintained online by Demaine, Mitchell, and O'Rourke HI.) 

Our Contributions. We describe a general proof technique that shows AAP-hardness of minimizing 
the stabbing number of perfect matchings, triangulations, and spanning trees. For the case of match- 
ings we show that it is also hard to approximate the minimum stabbing number within a factor below 
6/5. 

On the other hand, we present a mathematical programming framework for actually finding struc- 
tures with small stabbing number. Our approach characterizes solutions to stabbing problems as inte- 
ger programs (IPs) with an exponential number of cut constraints. We describe how the corresponding 
linear programming (LP) relaxations can be solved in polynomial time, providing empirically excel- 
lent lower bounds. Furthermore, we show that there always is an optimal fractional matching (or 
spanning tree) that contains an edge of weight above a lower bound of 1/3 (or 1/5 for spanning trees), 
allowing an iterated rounding scheme similar to the one developed by Jain for the generalized Steiner 
network problem ifTTl : compute a heuristic solution by solving a polynomial number of LPs. We 
conjecture that the objective function value of this heuristic solution is within a constant factor of the 
optimum. Our mathematical programming approach is also practically useful: as described in detail 
in our experimental study |[T2l . we can optimally solve stabbing problems for instances (taken from 
well-known benchmark sets of other geometric optimization problems) of vertex sets up to several 
hundred vertices. 

Our results in detail: 

• We prove that deciding whether a vertex set has a perfect matching of axis-parallel stabbing 
number 5 is an AAP-complete problem; we also extend this result to general stabbing number. 

'This should not be confused with the crossing number in graph drawing, which is the total number of crossing line 
segments. 
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• We prove that finding a triangulation of minimum axis-parallel stabbing number is an MV-haid 
problem; we also extend this result to general stabbing number. 

• We prove that finding a spanning tree of minimum axis-parallel stabbing number is an MV-haid 
problem; we extend this result to general stabbing number, and sketch AAP-hardness proofs for 
minimum axis-parallel or general crossing number. 

• We give an IP-based formulation for stabbing problems; the corresponding fractional LP solu- 
tions can be computed in polynomial time, providing a family of lower bounds. 

• We give results on the structure of fractional vertices of the resulting LP relaxation: for match- 
ing, we show that there always is an edge with weight at least 1/5, while for spanning trees, 
there always is an edge with weight greater than 1/3. This gives way to a heuristic algorithm 
based on iterated rounding; we conjecture that the resulting solution values are within a constant 
factor of the optimum. 

The vertex sets constructed in our hai^dness proofs make critical use of the coUinearity of vertices. 
On the other hand, our positive (LP-based) results do not make any assumptions on the structure of 
the vertex set: they can be used for vertex sets in degenerate as well as in general position, and can be 
applied to any family of stabbing hues that can be evaluated by considering a subset of polynomially 
many representatives. 

We have also performed a computational study on a diverse set of instances; the results show that 
our LP-based approach is good not only in theory (where we get a polynomial running time based on 
the ellipsoid method), but also for actually solving instances in practice (where we use the simplex 
method). Details are omitted from this theoretical paper; a report on the practical results can be found 
inHa. 

Related Work. Existing work dealing with structures of low stabbing number can be divided into 
algorithmic applications and implications on one hand, and extremal properties on the other hand. 

Agarwal fl] describes improved algorithmic solutions for problems such as ray shooting and im- 
plicit point locations queries (which by themselves have applications in polygon containment, implicit 
hidden surface removal, polygon placement, etc.); his main tool are spanning trees with low stabbing 
number. One of the theoretically best performing data structures for ray tracing in two dimensions is 
based on a triangulation of the polygonal scene; see Hershberger and Suri |rT6|: in their "pedestrian" 
approach to ray shooting, the complexity of a query is simply the number of triangles visited, i.e., cor- 
responds precisely to the stabbing number. Held, Klosowski, and Mitchell |[T5l investigate collision 
detection in a virtual reality environment, again, based on "pedestrian" ray shooting. More recently, 
Aronov et al. p| have performed an experimental study of the complexity of ray tracing algorithms 
and run-time predictors, which include average number of intersection points for a transversal line, 
and depth complexity. 

Extremal properties of crossing numbers were considered by Welzl |[26l and by Matousek ll20l . 
who showed that any planar set of n vertices has a spanning tree with a crossing number of 0{^/n), 
and provided examples requiring a crossing number of ^l{^/n). Another variant is studied by de 
Berg and van Kieveld [7|: the stabbing number of a decomposition of a rectilinear- polygon P into 
rectangles is the maximum number of rectangles intersected by any axis-parallel segment that lies 
completely inside of P; they prove that any simple rectilinear- polygon with n vertices admits a de- 
composition with stabbing number 0(log n), and they give an example of a simple rectilinear polygon 
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for which any decomposition has stabbing number Q{log n). They generaUze their resuhs to rectihn- 
ear polygons with rectihnear holes. Agarwal, Aronov, and Suri IH investigate extremal properties of 
the stabbing number of triangulations in three dimensions, where the stabbed objects are simplices; 
see also Aronov and Fortune i5j for this problem. Shewchuk |[24l shows that in d dimensions, a line 
can stab the interiors of ©(n^'^/^T) Delaunay d-simplices. This implies, in particular, that a Delaunay 
triangulation in the plane may have linear stabbing number. More recently, Toth ll25l showed that for 
any subdivision of d-dimensional Euclidean space, ti > 2, by n axis-aligned boxes, there is an axis- 
parallel line that stabs at least 0(log^/^'^~^^ n) boxes, which is the best possible lower bound. Gen- 
eralizations of the stabbing objects have also been considered: most notably, Chazelle and Welzl ||6l 
describe extremal properties of stabbing spanning trees in d-dimensional space by hyperplanes; the 
analogous problem and a corresponding result for matchings is also discussed in 161. 

This Paper. The rest of this paper is organized as follows. After some basic definitions and notation 
in Section m we give details of our various hardness proofs in Section [3] In Section HI we describe 
our LP-based approach for constructing bounds. Section [5]presents an iterated rounding technique for 
matching and spanning tree problems; we believe that the resulting algorithms yield constant-factor 
approximations. Final concluding thoughts and miscellaneous results and problems are presented in 
Section |6] 

2 Preliminaries 

Given a set L of line segments in the plane, the stabbing number of a line I is the number of segments 
of L that are intersected by I. The stabbing number of L is the maximum stabbing number over 
all lines t, the axis-parallel stabbing number of L is the maximum stabbing number over all axis- 
parallel lines I. In this paper, the set L will arise as the edges of a perfect matching, spanning tree, or 
triangulation of a given set P of n vertices in the plane, and our objective is to find such a structure of 
minimum stabbing number. Any reference to matching always means perfect matching. Therefore, 
when dealing with matchings, we assume that n be even, if necessary by omitting one of the vertices. 

We denote by St-Mat{P) the minimum stabbing number among all matchings of P, by St-Tre{P) 
the minimum stabbing number of all spanning trees of P, and by St-/\{P) the minimum stabbing 
number of all triangulations of P. We use St-Mat2{P), St-Tre2{P), and St-A2{P) for the minimum 
axis-parallel stabbing numbers. 

For a set L = {li, . . . , In} of line segments in the plane, the crossing number of a line £ is the 
number of connected components of UILi ^ ^- "^^^ crossing number of L is the maximum crossing 
number over all lines £; just like for the stabbing number, the axis-pai^allel crossing number is defined 
for axis-parallel lines. For matchings, trees, and triangulations, we use the analogous abbreviations 
Cr-Mat{P), Cr-Tre{P), and Cr-A(P), and their subscripted counterparts Cr-Mat2{P), Cr-Tre2{P), 
and Cr-/S.2{P) for the axis-paraUel crossing numbers. Note that stabbing and crossing number coin- 
cide for planar matchings. 

3 Complexity 

In this section we prove A/"'P-hardness of computing the minimum stabbing number of matchings 
and computing the minimum crossing number of triangulations; for spanning trees the proofs are 
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analogous, and we only give a sketch of the proof. Our technique is rather general and should be 
applicable to other structures and variants as well. 



3.1 Perfect Matchings 

Theorem 1 Deciding whether St-Mat2{P) < 5 is strongly MV -complete. 

Proof. Clearly, the problem is in MV. We show completeness using a reduction from 3SAT |[T3l . 
Assume we have a Boolean expression denoted by B{xq, x\, . . . , Xn-i) with n variables and k clauses 
of three literals each. We construct a set of vertices P that has a perfect matching M of stabbing 
number 5 if and only if the Boolean expression can be satisfied; in case of an unsatisfiable expression, 
a stabbing number of at least 6 cannot be avoided. 

Consider the overall layout of P as shown in Figure [T] We make critical use of the collinearity of 
vertices, using up all of the available stabbing number of 5 in a particular direction. Thus we are able 
to construct "barriers" which avoid any interference between the different gadgets. 
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Figure 1: Overall layout of the construction for St-Mat2{P). Shown is the layout for the 3SAT instance 
(xo V xi V xa) A (xo V X2 V X3) A {xi V X2 V X3), with a truth setting of xq = true, xi = true, X2 = 
false, X3 = true. Note that in full scale, the two gadgets representing the literal X3 in clauses and 1 
have disjoint x-coordinates. 



At the top of the layout are two groups of 10 vertices. The vertices in a group of 10 have the 
same y-coordinate. We call these two groups the top rows. The i-th vertex in the first top row has 
the same x-coordinate as the i-th vertex in the second top row. Below the top rows are 7i groups of 6 
vertices. All vertices in a group of 6 have the same x-coordinate, as shown in the figure. The vertical 
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lines through these groups of 6 vertices separate the variables from each other and from other vertices 
left of the variables. We call each such group a vertical barrier gadget. The vertices in banier i are 
sufficiently far below the vertices in baiTier i + 5 to ensure that horizontal lines through the vertical 
barriers have stabbing number at most 5. The barriers lie between vertical lines through the last and 
second last vertices of the top rows. To the left of the top rows and below the vertical barriers are 
k + 1 groups of 10 vertices. Each vertex in a group of 10 has the same y-coordinate. We call each 
such group a horizontal barrier gadget. The vertices in barrier i are sufficiently far to the right of 
the vertices in barrier i + 5 to ensure that vertical lines through the horizontal barriers have stabbing 
number at most 5. The horizontal barriers are used to sepai^ate clauses from each other, to separate the 
clauses from the variables, and to separate variables from other vertices above the variables. Between 
and to the right of the top two horizontal barriers are groups of 8 vertices. Each vertex in a group of 
8 has the same y-coordinate. We call each such group a reducer gadget of a variable. The vertices in 
reducer i are far enough to the right of the vertices in reducer i + 1 to ensure that vertical hnes through 
the reducers have stabbing number at most 5. 



true false false true 

column column column column 




Figure |2] shows a set of 6 vertices, numbered iq to in clockwise order, that represents a variable. 
These variable gadgets are separated by vertical barriers, and are placed below the vertical barriers, 
between the top two horizontal barriers and to the right of the reducers. Variables are placed such that 
the y-coordinate of the horizontal line through the reducer of a variable is in between the y-coordinate 
of vertices io and of that variable. The strip formed by vertical lines between vertices io and ii of 
variable Xi is called the Xi-column of the variable. The strip formed by vertical lines between vertices 
ii and 12 of variable xi is called the xi-column of the variable. 
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Figure 3: Clauses (xq V xi V X2) and (xq V xi V X2) with xq = X2 = true and xi = false. 
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The horizontal and vertical barriers create kn locations for literal gadgets. These ai^e groups of 
4 vertices representing the occurrence of a variable in a clause. Each group of 4 vertices forms an 
axis-pai^allel square. If a literal Xi appears in the clause cj, we place a literal gadget in the column 
of clause cj. If a literal Xi appears in the clause Cj, we place a literal gadget in the Xi column of clause 
Cj. The three literal gadgets for a clause are put on the same horizontal lines. The literal in column 
Xi of clause cj is to the left of the literal in column Xi of clause Ch for j < h. Similarly, the literal in 
column Xi of clause cj is to the left of the literal in column Xi of clause for j < h. Figure [3] shows 
the literals of two clauses. 

We first assume that B{xo,xi, . . . ,Xn-i) is satisfiable, and show that P has a matching M of 
stabbing number 5. We connect vertex i to vertex z + 1 in each of the two top rows for i = 0, 2, 4, 6, 8. 
We connect vertex i to vertex z + 1 in each vertical barrier for z = 0,2,4. We connect vertex i to 
vertex i + 1 in each horizontal barrier for i = 0,2, 4, 6, 8. We connect vertex i to vertex i + 1 in each 
reducer for i = 0,2, 4, 6. 

If the variable Xi has the value true, we connect the pairs (ic^s), {h,i2), and (^3,^4) of the 
variable. If the variable Xj has the value false, we connect the pairs {12, h)^ and {14,15). 

Notice that if Xi is true, then any vertical line in the -column stabs 2 edges in the top rows of M, 
and a vertical line in the Xj-column stabs 2 edges in the top row and 2 edges in the rectangle of 
the variable. If Xi is false, this situation is reversed. The column with vertical stabbing number 2 
is called the frwe-column of the variable, the column with vertical stabbing number 4 is called the 
false-column. In each literal gadget representing the value true we connect the four vertices with two 
horizontal edges. In each literal gadget representing the value false we connect the four vertices with 
two vertical edges. 

We can now verify that M has stabbing number 5. Any vertical line in the true-column of variable 
Xi stabs two edges in the top rows and at most two edges in a literal. Any vertical Une in the false 
column stabs two edges in the top rows, two edges of the variable and at most one in a literal. Any 
horizontal line in a clause stabs at most three literal gadgets, one of which is set to true. So these lines 
stab at most 5 edges of M. It can easily be verified that all other horizontal and vertical lines stab at 
most 5 edges from M. 

Conversely, we assume that P has a matching of stabbing number 5. We show that B is satisfiable. 
The matching used in this proof is illustrated in Figures [BE] and [3] Because the top rows contain 10 
vertices, these vertices have to be connected to each other, otherwise the stabbing number of P exceeds 
5. There are several ways to connect the sets of 10 vertices. If we connect vertex i to vertex i + 1 
for i = 0, 2, 4, 6, 8 in each row, the number of edges stabbed by any horizontal or vertical line is 
minimized. Therefore we may assume without loss of generality that these edges are in the matching 
M. Collinear vertices are dealt with in a similar manner: If there are several ways to connect a set 
of collinear vertices, we will prefer connections that have minimal stabbing number for all stabbing 
lines. 

Thus, we can connect vertex i to vertex i + 1 for z = 0, 2, 4, 6, 8 in each horizontal barrier gadget. 
For the same reason, we can connect vertex i to vertex z + 1 for i = 0, 2, 4, 6 in each reducer gadget. 
Because vertical lines through the vertical baiTier gadgets stab two edges in the top rows, we can 
connect vertex i to vertex i + 1 for z = 0, 2, 4 in each vertical barrier gadget. 

Now it is easy to see that no vertex in a variable or literal gadget for some variable Xj can be 
matched with any vertex not involved with representing the same variable Xi'. otherwise, the edge 
would cross a line thr^ough a vertical barrier, which akeady crosses another five edges. Furthermore, 
such a vertex must be matched with vertices from the same gadget (either the same variable gadget or 
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the same literal gadget): otherwise, we get a violation at a line through a horizontal banier gadget. 

Now each reducer gadget contributes four to a horizontal stabbing number. Thus, we cannot 
connect the six vertices of a variable Xj by three vertical edges. Figure |2] shows the two remaining, 
essentially distinct matchings. 

Of the three literals in each clause, one has to be set to true, otherwise there will be a horizontal 
stabber intersecting six edges. The true literal, say, x, must lie in a true-column of a variable, because 
vertical lines in this column have stabbing number two. Any other literal in this column can also be 
set to true. The literals x lie in the/a/^e-column of the same variable, and have to be set to false. So 
if a matching of stabbing number five exists, there is a truth assignment of the Boolean expression. □ 

Corollary 2 There is no a-approximation algorithm for St-Mat2{P) with a < 6/5; in particular 
there is no polynomial time approximation scheme (PTAS), unless V = NV. 

Corollary 3 Computing St-Mat{P) is a weakly NV-hard problem. 

Proof. We apply a perturbation technique, similar to the one in lITOl . We start with the same basic 
construction as for the hardness proof for the axis-parallel case, and consider the grid formed by the 
coordinates of the resulting vertex set. This grid is modified such that the interpoint distances between 
the vertices of the same gadget are B(e" ~'^^) for the literal gadgets and 0(e" ) for all other gadgets. 
Furthermore, the rest of the grid is perturbed by powers of e, such that only axis-parallel lines can stab 
more than two gadgets; in particular, we increase the vertical distance between variable gadgets and 
the (narrower) literal gadgets by a sufficient amount, in order to make sure that no line through two 
literal gadgets for the same variable can intersect the corresponding variable gadget. Now it is easy to 
see that lines that are not axis-parallel can stab at most four line segments, leaving only axis-parallel 
lines as critical. □ 

3.2 Triangulations 

Our basic proof technique is the same as for matchings. We first describe the construction of barrier 
gadgets, using the following terminology. A horizontal line is given by a set of vertices that are 
horizontally coUinear. A vertical line is given by a set of vertically coUinear vertices. A row consists 
of two horizontal lines, and the (empty) space between them. A column consists of two vertical lines, 
and the (empty) space between them. 

Lemma 4 Consider a row consisting of two horizontal lines la and If, in P, having a and b vertices, 
respectively. If the combined number of edges on la and lb is a + b — i — 2, then a horizontal stabber 
between la and lb encounters at least a + b + i — 2 triangles in any triangulation of P and its crossing 
number is at least a + b + i — \. 

Proof. Assume without loss of generality that la lies above lb- Suppose there are a, — ia — 1 edges on 
the line la and b — % — 1 edges on the line lb with ia + ib = i- For each edge {u, v) on there is a 
triangle (u, v, w) where w lies either on or below lb. Let A denote this set of triangles. Similarly, for 
each edge {u, v) on lb there is a triangle {u, v, w) where w lies either on or above la. Let B denote 
this set of triangles. For each two neighboring vertices u and v on la for which there is no edge (n, v) 
there ai^e triangles (u, uo,ui) and {v, vo,vi) such that uq and vq lie on or below lb and ui and vi lie 
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above la- Let la denote this set of triangles. Also for each two neighboring vertices u and v on for 
which there is no edge (u, v) there are triangles {u, uo,ui) and {v, fo, ui) such that uq and vq lie on 
or above and ui and vi lie below li,. Let denote this set of triangles. It is not hard to verify that 
any two of the the four sets of triangles A, B, la and If, have an empty intersection. A horizontal line 
I between la and /f, stabs every triangle in A, B, la and I^,. So / stabs at least \A\ + \B\ + \Ia\ + \Ib\ 
= {a — ia — 1) + {b — if) — 1) + 2ia + 2if, = a + b + i — 2 triangles and a + b + i — 1 edges. □ 

The lemma holds analogously for two vertical lines that form a column. When a row consists of 
two horizontal lines that have Cr-/S.2{P) + 1 vertices altogether, we call it full or fully triangulated. It 
follows from the lemma that all Cr- A2 (-P) — 1 edges on the Unes and lb have to be present. 

Theorem 5 Finding Cr-A2{P) is MV -hard. 

Proof. Again we use a reduction from 3SAT, and the proof proceeds along the lines of the proof of 
Theorem [U See Figured for a schematic layout of a representing vertex set P for the 3SAT instance 
B{x 0, xi,X2) = {xq V xi V X2) A (xo V xi V X2) A (xq V xi V X2)- Figure [5] shows the structure of 
variable gadgets. 




Figure 4: Overall layout for Cr-A2{P). Clauses are (xq Vxi VX2), (xq Vxi VX2), and (xq Vxi VX2), 
with Xq = false and xi = X2 = true. An^ows indicate full rows and columns, light or dark shading 
indicates true or false variables and literals as before. 
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For a given Boolean expression B{xq, xi, . . . , with n vaiiables and k clauses of three lit- 

erals each we construct a set P of vertices. We show that there is a value K such that B is satisfiable 
only if Cr- A2(P) =2K -l\ iiB cannot be satisfied, Cr-A2(P) is at least 2K. 

In Figure |4] we have K = 39 and a grid of vertices with some well-defined holes. The maximum 
number of vertices in a horizontal or vertical line is K, and many lines have exactly K vertices. By 
Lemma m a full row or column in this setting has exactly 2K — 2 triangles. Gadgets are separated by 
full rows and columns. 



false true 
column column 




K -I 



Figure 5: A variable gadget and how it is embedded in a grid of vertices. Arrows indicate full rows. 



Figure [5] shows two horizontally aligned rectangles of eight vertices each that together represent 
a variable Xj. We call the strip formed by vertical lines that stab the left rectangle the Xi-column, 
and the strip formed by vertical lines that stab the right rectangle the Xi-column of the variable. The 
gadget works essentially the same way as that in Figure |2] Each rectangle has full rows and columns 
as neighbors. We indicate how this can be achieved in horizontal direction in Figured By Lemma|4] 
we conclude that all edges on the boundary of the convex hull of each rectangle ai^e present in any 
triangulation of minimal crossing number. The horizontal lines that contain the top and bottom lines 
of the two rectangles, respectively, contain K vertices each; the horizontal line that passes through 
the middle of the rectangles contains K — 1 vertices. Therefore Lemma H] shows that exactly one 
edge along this horizontal line in the middle of the rectangle may be missing. We call the strip that is 
spanned by all vertical lines that stab the rectangle with the missing horizontal edge the true-column 
of the variable Xj. The strip that is spanned by all vertical lines that stab the rectangle for which the 
middle horizontal edge is present is the /aZ^e-column of the variable Xj. As shown in Figure IH any 
satisfying literal adds exactly one less to the vertical crossing number than an unsatisfied one. In the 
overall layout, variable Xj is placed below and to the left of variable Xj for i < j in such a way that 
variables ai^e vertically separated from one another by a full row, and horizontally separated from one 
another by a full column. 

Each literal is represented by a square with eight vertices on its boundary. We make the width of 
rectangles of the variables equal to a power of two, and wide enough to accommodate the necessary 
number of literals. That is, each rectangle of a variable is of a width at most four times the number of 
occun^ences of the most frequent literal in B. Figure [5]gives a hint at how this widening of a rectangle 
is done. Notice that if the most frequent literal in B occurs t times, there are ©(logt) rows above 
and below the variable so that in the top and bottom lines of these groups of rows the vertices above 
and below the variables are a distance of one apart. The three literals of a clause Cj are horizontally 
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aligned, and the two rows that are spanned by them ai^e called the clause cj. Clauses ai^e sepai^ated 
from each other by full rows. If a literal Xi appears in the clause Cj, we place a literal gadget in the 
Xj-column of the clause gadget Cj. If a literal Xi appears in the clause Cj, we place a literal gadget in 
the X j-column of the clause gadget Cj. The literal in column x i of clause Cj is to the left of the literal 
in column x j of clause for j < h. Similarly, the literal in column x j of clause cj is to the left of 
the literal in column X j of clause for j < h. So no vertical line stabs the interior of more than one 
literal. 

By adding vertices to the right of the literals we ensure that a horizontal line through the top or 
bottom row of the three literals of a clause has exactly K vertices, and a horizontal line through the 
middle horizontal line of the three literals of a clause has exactly K — 2 vertices. So two edges along 
these middle lines may be missing in a triangulation of minimum crossing number, but no more than 
two. As we will argue later, these edges will be missing in the interior of at most two of the literals. 
We call the missing of the horizontal middle edge in a literal the false setting of the literal, and the 
presence of this edge within a literal the true setting of the literal. 

Because the rows above and below a variable are full, we can assume that they are triangulated 
as shown in Figure \5\ because any other triangulation would result in strictly higher vertical crossing 
numbers. So the true column of a variable has a vertical crossing number that is one less than the 
crossing number of a. false column. In adding vertices at the bottom of the clauses we ensure the 
following vertical vertex counts. First of all, the columns neighboring the variables have to be full; in 
particular, vertical lines that stab the left or the right vertices of a variable rectangle have K vertices 
each. All other vertical lines through a variable should get a number of vertices so that if the corre- 
sponding column of the variable is set to false, all remaining edges on these vertical lines have to be 
present. This implies that if the column is set to true, we can have one missing edge in the vertical 
line to the right of the left boundary of the rectangle, and if this edge is missing, there cannot be a 
missing edge in the next vertical line, one edge missing in the next line, etc. 

Let B{xo,xi, . . . , Xn-i) be satisfiable. We show that P has a triangulation of crossing number 
2K — 1 that is minimum by LemmalU All full rows and full columns are fully triangulated. If variable 
Xi has the value true, we triangulate the interior of the two rectangles of variable Xj in such a way 
that the Xj-column becomes this variable's true-column, and the Xj-column becomes this variable's 
false-column. The triangulation of the interior of the rectangles is reversed when variable Xi has the 
value false. We set each literal that represents the value true to its true setting, and set each literal that 
represents the walue false to its false setting. The triangulation can be completed arbitrarily. 

We can now convince ourselves that such a triangulation of P indeed attains Cr-A2{P) = 2K — 1. 
No fully triangulated row or column has a crossing number larger than 2K — 1. Because exactly one 
edge is missing in the horizontal middle line of each variable. Lemma |4] implies that the crossing 
numbers of the two rows of a variable are both equal to {2K — 1) + 1 — 1 = 2K — 1. In each clause 
there is at least one literal in its true setting. Therefore, we can afford two extra triangles caused by the 
false setting of the other two literals in the clause, and no row of a clause intersects more than 2K — 2 
triangles. Finally, we may have one edge missing on the vertical lines passing trough the middle of a 
literal only in a true-column. This condition holds by definition of the setting of the literals according 
to the truth value of the variables. The "arbitrary completion" of the triangulation only happens in the 
lower right of the construction. In this corner, vertical and horizontal lines have a low vertex count 
(except for the boundary), and the allowed crossing number is not exceeded. 

For seeing the converse, assume that there is a triangulation of P that has crossing number 2K — 1. 
We show that B is satisfiable. Because a full row or column can be triangulated in such a way that 



11 



the crossing number of 2K — 1 is not exceeded, we only have to take care of the rows and columns in 
which we have a degree of freedom, and where the vertex count is critical by Lemma|4] Because both 
the top row and the bottom row of a clause have 2K — 2 vertices, we can afford at most two literals 
of each clause set to false. One literal in each clause has to be set to true. The true literal has to be in 
a true-co\\mm of a variable, for otherwise the vertical crossing number would exceed 2K — 1. Any 
other literal in this column can also be set to true. The horizontal vertex count of the horizontal lines 
of the variable forces the second column of this variable to be a/a/^e-column. In order not to exceed 
the allowed crossing number in vertical direction, all literals in this column have to be set to false. As 
one easily checks, this yields a consistent setting of the variables, and B is satisfiable. 

Finally, the size of our construction is indeed polynomial: let n and c be the number of variables 
and clauses of B{xo, xi, . . . , Xn-i). Let t be the number of times the most frequent literal in used in 
B. A rectangle that represents a variable has width at most 4t and requires ©(log t) rows. The number 
of rows used by the clauses is 0(c). In order to achieve the correct vertex count in each line we may 
have to add vertices to the right and below, which requires at most B(t) additional rows and columns. 
Therefore K is polynomial in n, c and t. □ 



3.3 Spanning Trees 

The basic construction for showing hardness of finding a spanning tree of minimum stabbing number 
is similar to the one for matchings. As before, we use barriers to restrict possible connections: we 
make use of the arrangement shown in Figure [6l which works as a bamer gadget because of the 
following lemma. 



Lemma 6 Consider three parallel lines, l\, ^2. ^3. with a set Si ofk vertices on line ii, i = 1,2, 3; 
let S = S*! U §2 U 53. Consider P D S and a spanning tree T of P with stabbing number k+\. Then 
no edge ofT crosses the strip spanned by the three lines. 




■\ k vertices 



(a) 



(b) 



Figure 6: A horizontal barrier gadget, consisting of 3k vertices: (a) In a spanning tree of stabbing 
number + 1, no edge may cross the shaded region, (b) Symbol for the baiTier gadget; the dotted line 
indicates the blocked strip. 

Proof. See Figure [6l Consider a spanning tree T of P, with w G P \ 5 lying outside of the strip. 
Orient all edges of T towai^ds v. Each vertex in S must have outdegree 1, meaning that there are k 
outgoing edges for each of 51,52,53, contributing k to the stabbing numbers along £1, ^2, ^3- One of 
the outgoing edges of 52 must intersect ii or ^3 in order to connect 52 to the rest of the graph; thus, 
one of those two lines stabs k + 1 edges, implying the claim. □ 

Our variable gadgets look as in Figure |71 shown is the gadget for variable x,; note that the 
gadgets for xq,. . . ,Xi-i are left and below the box spanned by the gadget, while the gadgets for 
Xj+i, . . . , Xn-i are above and to the right of the spanning box. The bold squares below the gadget 
indicate the position of literal gadgets, which wil be discussed further down. 
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Also note the use of vertical barrier gadgets: a number of i next to it indicates a gadget consisting 
of 3i vertices, which akead requires a crossing number of i + 1; thus, only k — i additional edges in a 
spanning tree may cross the dotted line induced by such a gadget. The aiTows pointing down from the 
bottom indicate a number of literal gadgets, consisting of 2 x 2 arrangements of vertices. See Figured] 
for the resulting overall arrangement. 

Lemma 7 Let S be the arrangement of vertices shown in Figure^ with barrier gadgets placed and 
sized as indicated, and let P ^ S. Let P be constructed as shown in Figure |S] Then any spanning 
tree of P that has stabbing number at most k + 1 must use at least one of the two edges at the bottom 
of the arrangement, labeled ej'^ (for truej or el^ (for false.) 




Figure 7: A variable gadget for vaiiable Xj. (a) In a spanning tree of stabbing number k + 1, the true 
or the false setting is chosen, (b) Symbol for the variable gadget. 

Proof. Assume there is a spanning tree of stabbing number at most k+1. Consider the barrier gadgets 
labeled Bl, B2, B\, B\, and the corresponding lines, £[, C^^ i\, £\. By the previous lemma, no edge 
can cross one of those lines. Therefore, the literal boxes below each clause must be connected within 
the vertical strips bounded by £2 and £| on one hand, or £| and i\, on the other hand. This requires at 
least one edge within each of the two strips to cross the line £q. Moreover, the lines £\ and i\ must not 
be crossed, implying that the variable gadgets are connected to their neighbors at barriers B\ and B^. 
These connections between neighboring variable gadgets form a stair-like chain of variable gadgets, 
as shown in Figure (8) connections ai^e in increasing x- and y-order, which con^espond to increasing 
variable indices; we call this the exterior stair. Similarly, the barrier gadgets associated with each of 
the two truth settings of a variable form one stair-like chain each, as shown in Figure |7J we call these 
the two interior stairs. 
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Now consider the horizontal barrier Bq, consisting of three groups of A: — 5 — 2(n — 1 — i) vertices 
each. By the previous arguments, the line £q has to cross all of the edges connecting the true and false 
literal boxes of the n — l — i variables with higher indices, i.e., cross 2(n — edges. Furthermore 
in variable xi there are four other edges that connect vertices above line /q to vertices below /q and so 
are crossed by /q. For example these edges could be the ones labeled c^, Cg, c\,c\ in the figure. This 
allows only one of the edges ej'^ and ej^ to be used for connecting the two interior stairs with each 

other and with the exterior stair; thus, at least one of the edges ej^ and ej'^ must be used, proving the 
claim. □ 

Making use of the above gadgets, we get 
Theorem 8 It is MV-hard to determine St-Tre2{P)- 



Proof. The basic idea for the construction is similar to the one used in the previous sections, making 
use of Lemmas [6] and |7] The use of gadgets and the overall layout of the construction are shown in 
Figure [U 




Figure 8: The overall layout for the hardness proof for spanning trees. There is a total of n variables; 

is a sufficiently large number. In question is the existence of a spanning tree with stabbing number 
fc+1. Shown is the representation of the 3SAT instance {xqVxi VX2) A (xq Vxi VX2) A (xq Vxi VX2), 
for n = 3, with xq = false and xi = X2 = true. 

Given a Boolean expression B{xq, xi, . . . , Xn-i), we can find a vertex set that has a spanning 
tree of a stabbing number A: + 1, if and only if B is satisfiable. Consider the vertex set as given in 
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Figure [8] with k = 3?i. If B{xq,xi, . . . ,Xn-i) has a satisfying truth assignment, we first connect 
the vertex on the left side of the drawing into one long path. This path contains all horizontal and 
vertical bamer gadgets as shown in Figure [8] We connect the variable and literal gadgets according to 
their values in B. In each true- and/a/^e-column, there is a vertex lower than all horizontal barriers, 
and to the right of all literals in that column. We connect this vertex to the right lowest vertex of the 
lowest literal in the same column. Because each clause has at least one true literal, the horizontal 
stabbing number of a stabber through a clause is at most k, because the horizontal barrier on the left 
contributes k — 2n — 1, the variable(s) set to true and the n — 3 variables that do not appear in the 
clause contribute two each, and the vaiiable(s) set to false contribute thi^ee each, which is at most 
k — 2n — 1 + 2{n — 2)+6 = A;+l. A horizontal stabber through the horizontal bamer in the i-th 
variable gadget stabs k — 5 — 2{n — 1 — i) + 1 edges in the barrier, one more edge to the left of the 
barrier, four edges of the variable to the right of the gadget and two more for each of the n — i variables 
to its right, for a total of A; — 5 — 2(n — 1 — z) + 1 + 5 + 2(n — 1 — i) = k + l. A vertical stabber 
through literals stabs at most k — 1 edges in the vertical baniers, plus two more, either one from a 
false literal and one in the variable gadget, or two edges in the true literals. So we have a spanning 
tree of stabbing number A; + 1. 

Conversely, assume that the vertex set has a spanning tree of stabbing number k + 1. As we 
showed in the proof of Lemma |7J the literal gadgets can only be connected within their respective 
strips, forcing at least 2?i stabbed edges. Furthermore, each false literal (being connected in a "u"-like 
fashion) causes an additional stabbed edge, while a true literal (connected in a "c"-like fashion) does 
not cause any additional stabbings. Thus, each clause must contain at least one true literal. Because 
of Lemma U\ at least one of the edges e j. and e} must be present, guaranteeing that only the negated 
or only the unnegated literals for each variable can be connected in a "c"-like fashion, i.e., forcing a 
feasible setting of the variables. Thus, we get a truth setting of the variables that satisfies B. □ 

This immediately implies 
Corollary 9 It is MV-hard to determine St-Tre{P). 

Proof. The argument is similar to the one in Corollary [3] for matchings: use the construction of 
Theorem m for which the criticality of certain axis-parallel lines requires satisfying a 3SAT instance 
in order to achieve low stabbing number. Scale down the bounding boxes for all gadgets, with literal 
gadgets ending up in appropriately smaller bounding boxes. Then perturb the position of the gadgets, 
shifting all vertices in the same gadget by the same amount, such that no line can intersect the bound- 
ing boxes of any three gadgets, again inserting a sufficient vertical distance for excluding a line that 
stabs a variable gadget and two of its literal gadgets. This leaves only the axis-parallel lines to be 
critical, implying the same combinatorial behavior as in the axis-parallel case. □ 

The hardness proof for minimizing the crossing number has the same structure as the one for 
stabbing number. Instead of the barrier gadget implied by Lemma [6l we use a slightly different one, 
as shown in Figure |9l 

Lemma 10 Let S be the kx((k — 1)^ + k) arrangement of vertices shown in Figure^ and let P ^ S. 
If P has a spanning tree T with crossing number k, then no edge ofT connecting two vertices outside 
of the arrangement crosses the horizontal strip spanned by the arrangement. 

Proof. Suppose there is a crossing edge. Consider the (fc — 1)^ + A; horizontal lines that pass through 
k vertices of the vertex arrangement. Because of the crossing edge, the intersection of each line with 
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Figure 9: A barrier gadget for showing hardness of minimizing the crossing number of a spanning 
tree. 

the edges and vertices within each corresponding horizontal set cannot consist of more than k — 1 
connected components, so this intersection must contain at least one horizontal edge within the set; 
thus, a total of at least {k — 1)^ + A; horizontal edges within the arrangement cannot be avoided. Now 
consider the k vertical lines that pass thi^ough {k — l)'^ + k vertices of the aiTangement; these lines 
subdivide the plane into k+l vertical strips, k—l of which have width 1. Any of the at least (fc— 1)^+A; 
horizontal edges within the arrangement must cross at least one of the intermediate strips; thus, the 
average number of horizontal edges per bounded strip is at least ^^"^"'j-^"^^" = (fc — 1) + > k. 
By the pigeonhole principle, this implies that there must be a bounded vertical strip that is crossed by 
more than k edges, a contradiction to our assumption that the crossing number is at most k. □ 

As the figure shows, there is a feasible subtree of the gadget, as long as no edge crosses the 
indicated strip. This allows us to use the aiTangement as a barrier gadget. As for the rest of the 
construction for the proof of Theorem [8l stabbing and crossing numbers coincide, this immediately 
implies the following. 

Theorem 11 It is MV-hard to determine Cr-Tre2{P). 

From this it is easy to derive the following, again using a perturbation argument. 
Corollary 12 It is MV-hard to determine Cr-Tre{P). 

Proof. The argument is similar to the one in Corollary |9l use perturbation to ensure that only axis- 
parallel lines can stab more than two gadgets. The only additional difficulty that has to be overcome 
is the fact that a diagonal line through one of the barrier gadgets in Figure |9] may have crossing 
number k: in principle, crossing number 2k could arise from stabbing two such gadgets in a non-axis- 
parallel manner. However, scaling the bounding boxes of the bamer gadgets in a way that the vertical 
dimension is much smaller than the horizontal dimension (say, by a factor of 0(e" )) makes sure that 
only Mnes of slope within range [—e^^,e^^] achieve crossing number k for one gadget. Thus, only 
almost horizonal lines are of concern; it is easy to see that the above perturbation guarantees that no 
such line intersects more than one barrier gadget. Therefore, only axis-parallel lines can be critical, 
and the hardness proof remains valid. 
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4 Integer Linear Programs for Minimum Stabbing Number 



In view of the negative complexity results for our problems there are two major directions to proceed: 
providing (good) lower bounds on the minimum stabbing number in order to obtain approximation 
algorithms; and insisting on optimality despite AAP-hardness. Our (integer) linear- programming ap- 
proach is an elegant way to combine both issues. We deal with them in the next two sections. 



4.1 Perfect Matchings 

In combinatorial optimization, P corresponds to the vertex set of a straight-line embedded complete 
graph G = {P, E); then a matching M can be represented by its edge incidence vector x G {0, 1}^, 
where xij = 1 if ij G M, and Xij = otherwise. Using these variables we are able to state an integer 
lineal- program for finding a perfect matching of minimum stabbing number. For S C P, denote by 
5{S) = {ij G E \ i £ S,j ^ S} the cut induced by S. 

minimize k (1) 

s.t. Xij = 1 Vi e P (2) 

ijesai}) 

Xij>l V0 / 5 c P, |5| odd (3) 

Xij < k V stabbing line £{d) in direction d (4) 

Xij €{0,1} yijeE (5) 
We obtain the associated linear programming (LP) relaxation by replacing ^ by 

Xij > dll) 

The inequalities (O and (O ai^e necessarily satisfied for any perfect matching given by x, where the 
blossom constraints (O ensure that any subset of vertices of odd cardinality has at least one edge to 
the outside. In his seminal paper Edmonds Q showed that — together with — these inequalities 
already constitute the complete description of the perfect matching polytope, that is, its extreme ver- 
tices exactly correspond to the incidence vectors of perfect matchings in G, i.e., the LP relaxation is 
integral. 

To this description we add the stabbing constraints for which we count the number of inter- 
sections of matching edges with any given line; this number is bounded by the variable k, and k is 
minimized. We have to choose this way of modeling because of our min-max objective. An opti- 
mal solution X to the integer program ([T|)-(l5l) represents a matching with stabbing number exactly 
St-Mat{P). In a pure integer programming description, the blossom constraints ai^e implied; however. 
Figure [TO] shows that the introduction of stabbing constraints yields a polytope that is no longer inte- 
gral, which is to be expected for an MV-haid problem. That is why we consider the LP relaxation, 
for which the use of blossom constraints yields a considerably tighter set of solutions. We will make 
use of this fact in the following section; here we comment on the complexity of the given LP and its 
solution. 

We can limit the number of stabbing constraints (0]): for the axis-parallel stabbing number, we 
can assume without loss of generality that a critical stabbing line runs through a vertex, so we get 
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Figure 10: An optimal fractional solution of value 4/3 with maximum edge weight 1/3. 



at most 0{n) restrictions for the axis-parallel stabbing number; for the general stabbing number, we 
may assume that a critical stabbing line runs through two vertices, so we can limit ourselves to O(n^) 
restrictions. On the other hand, we have exponentially many blossom inequalities; however, it is 
well known that blossom inequalities can be separated in polynomial time |[22l . i.e., one can check 
in polynomial time whether a given x violates some blossom inequality, and if so, identify such an 
inequality. This polynomial-time separation allows us to solve the linear programming relaxation ([T])- 
dsTI) in strongly polynomial time |[23l Thm. 5.11] by means of the ellipsoid method |[T4l Thm. 66.5]. 
An optimal solution x will in general be fractional, and we speak of fractional stabbing number in 
this context. It is a lower bound on St-Mat{P). 

4.2 Spanning Trees 

There are several polynomial-size LP formulations for spanning trees, see e.g., |fT9l|. However, similar 
to matchings, we choose an exponential-size integer program that is again based on cut constraints. 
We directly state the LP relaxation. 

minimize k (6) 

s.t. ^ Xij = n — 1 (7) 

ijeE 

^ Xij>l V0 / 5 C P (8) 

ij£5{S) 

Xij < k V stabbing hne £{d) in direction d (9) 

ij:ijne{d)^i 

X > (10) 

Equation dT]) ensures the right number of edges in a tree solution, and connectivity is given by the cut 
constraints Just like blossom constraints for matching problems, separation over these constraints 
can be done in polynomial time by means of a minimum cut routine, so this LP can also be solved in 
strongly polynomial time. 
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5 Iterated Rounding 



The above LP relaxations provide lower bounds for optimal solutions, but no upper bounds in the 
form of feasible solutions for our stabbing problems (except if we solve the integer programs, which 
may take exponential time). Therefore, our next objective is to find, in polynomial time, an integer 
solution that is not too far from an optimal one. 

For this purpose we consider the support graph of a (fractional) solution x: it consists of all edges 
e that have a strictly positive value Xe > 0. If we could be assured that in any optimal LP solution, 
all edges e in the support graph solution were 6-heavy, i.e., x^ > h for some constant & > 0, then 
we could consider rounding up all positive edge weights to one in order to get a l/&-approximation. 
Unfortunately, Figure [TT] gives an indication that there may not be such a lower bound. 

An alternative to simple LP rounding is offered by the iterated rounding technique introduced by 
Jain ifTTl for generalized Steiner network problems: at each step, it suffices to identify one heavy edge, 
round it up to one, fix it, and re-solve the remaining LP. The key ingredient is that the modified LP 
also has a heavy edge, so the process can be iterated, hence the name: iterated rounding. 

A crucial step in Jain's seminal work was the fact that for so-called generalized Steiner networks, 
1/2 is a valid lower bound for the value of the heaviest edge of fractional weight. A sufficient condition 
for the validity of this bound that is purely combinatorial is that all constraints in the LP relaxation are 
cut constraints, i.e., of type Ylie(^5{s) — ^{^)- 

Unfortunately, stabbing constraints do not fall into the category of cut constraints, and it can be 
seen from Figure [TO] that Jain's lower bound does not hold for stabbing problems. This makes it 
necessary to establish a separate lower bound on the value of the heaviest fractional edge. We will do 
this by exploiting the underlying geometric nature of our optimization problem, and establish the fact 
that there always is an optimal LP solution with planar support graph. Planarity is proven by shifting 
weight from a crossing pair of edges to a non-crossing one; for matchings (Lemma [T3l) this requires 
some extra care because of the blossom inequalities. The proof for spanning trees (Lemma O is 
almost completely analogous. 

Lemma 13 For any even set of vertices in the plane, there is a fractional perfect matching x of 
minimum stabbing number, such that the support graph ofx is planar Such a fractional matching can 
be found in polynomial time. 

Proof. The set of all LP solutions is bounded and nonempty, so the set of all optimal solutions is 
nonempty and compact. From this set, consider a solution x that minimizes the total stabbing number, 
i.e., the sum of stabbing numbers over all combinatorially different lines. We claim that the support 
graph of x cannot contain any crossing pair of edges. Refer to Figure [121 

Suppose there were ei3 := {fi,f3}, 624 := {u2i^^4} with Xej3 > and Xe24 > 0, such that 613 
and 624 cross. We will argue that this implies the existence of an alternative solution of the same 
maximum stabbing number, but infinitesimally smaller total stabbing number, i.e., a contradiction to 
our assumptions on optimality. 

Consider ei2 := {^1,^2}, 634 := {f3,f4}, £14 := {fi,f4}, 623 := {^2,^3}- Among all LP 
constraints (fB-dsTl). let s{x) be the smallest positive slack (i.e., difference between left-hand and 
right-hand side), and choose < e < s(a;). As Ylie&5{vi) = 1 and Xe^^ > £, Xe24 > s, we have 

Consider two possible alternative solutions arising from shifting an e of weight from e 13 and 624. 
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Figure 11: Finding a fractional matching of minimum stabbing number: shown is an optimal LP 
solution for an instance of 100 random vertices; see |[T2l for technical details and results for other 
instances. 
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Figure 12: Improving a fractional solution with crossing edges while preserving all blossom inequal- 
ities. 



Let x' and x" be defined by 




e e {ei3, 624}, 

eG {612,634} (11) 
all other e. 



and 




6 G {613,624}, 

Xg := ■{ Xe + e for 6 €{614,623}, (12) 

all other e. 

By convexity, both x' and x" satisfy all stabbing constraints that are valid for x. Furthermore, it is 
easy to see that both x' and x" have smaller total stabbing number than x. Thus, it suffices to argue 
that at least one of these solutions satisfies all blossom inequalities. 

Assume that for each of the two alternative solutions, a blossom inequality is violated; by our as- 
sumption on £, this means that there are two odd sets. Si d P and ^2 C P, such that X^gg^^^^) x'{e) < 

1 and Eee5(52) ^"(e) < 1' ^hile Eee5(Si) = 1 and EeG5(52) = 1- Let Si := \{vi,V2,V3,V4}n 
Si\. It is straightforwai^d to see that for s j G {0, 1, 3, 4}, both x' and x" satisfy all blossom inequalities 
that are valid for x; similarly, the blossom inequalities are satisfied for v 1, ^3 G ^i and ^2,^4 G Si, as 
See5(Si) ^'{^) > Z]eG5{Si) ^{^)- Therefore, we only need to consider v i,V2 G 5i and ^3,^4 G Si, 
and vi,Vi G ^2 and V2,V3 G Let Ti := Si n ^2, T2 := SiH^, :=^n'S^T4 :='S^nS2, 
and Eij ■= {e = {x,y} e E \ x e Ti,y e Tj}. As \Ti U Ta], |Ti U Ti], \T2 U r3|, and IT3 U T4I 
are all odd, we may assume without loss of generality that |Ti| and IT3I are odd, and IT2I and {T^l are 
even. 
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Then we have 



2 = ^ xe+ ^ Xe (13) 

ee5{Si) ee5(S2) 

= ^ Xe+ ^ Xe + 2 ^ Xe (14) 
eG5(Ti) ee<5(T3) e(^E24 

- ^ Xe + ^ Xe + 2Xe24 (15) 
ee<5(Ti) eG<5(T3) 

> ^ Xe+ ^ Xe- (16) 

ee<5(Ti) ee<5(T3) 

This implies that minj^gg^^j,^^ Xg, X]e6(5(r3) ^e} < 1> contradicting the assumption that x satisfies 

all cut inequalities. □ 

The above lemma establishes the existence of optimal non-crossing matchings; when actually 
trying to find a non-crossing matching, it suffices to consider an additional term in the objective 
function of our LP that refers to the total length of the edges. This ensures that the total length 
of matching edges is minimized, avoiding crossings in the first place. For more details, see our 
experimental paper |[T2l . 

For the case of spanning trees, we can use a similar approach, based on a similar lemma. 



Lemma 14 For any set of vertices in the plane, there is a fractional spanning tree x of minimum 
stabbing number, such that the support graph of x is planar. Such a fractional spanning tree can be 
found in polynomial time. 



Proof. We proceed completely analogous to the proof of the previous lemma to deduce that if there 
was a fractional solution x of minimum total stabbing number among all solutions with optimal stab- 
bing number, such that the support graph of x has two crossing edges, then we could shift an in- 
finitesimal amount of weight from the crossing edges to two non-crossing ones, such that the stabbing 
number remains the same, but the total stabbing number decreases. 

The technical steps of the argument are virtually identical: consider two crossing edges ei3 and 
624 as shown in Figure [l2l and the alternative solutions Xg and Xg as defined in (fTT]) and (fT2l) . For 
arguing that at least one of these solutions is feasible, we do not have to consider blossom constraints, 
but cut constraints of type (HJ). 

Assume that for each of the two alternative solutions, a cut constraint is violated; by our assump- 
tion on e, this means that there are two sets, C P and 52 C P, such that x'(e) < 1 and 

Eee5(S2) ^"(e) < 1' while EeG<5(5i) = 1 and EeG5(52) = 1- Let Si := \{vi,V2,V3,Vi} n Si\. 

It is straightforwai^d to see that for Sj G {0, 1, 3, 4}, both x' and x" satisfy all cut inequalities that 
are valid for x; similarly, the cut inequalities are satisfied for u 1,^3 G 5*1 and f2,f4 G ^i, as 
SeG(5(Si) ^'{^) > J2e£5{Si) ^(^)- Therefore, we only need to consider v i,V2 G -Si and vs,V4 G Si, 
and vi,Vi G S2 and V2,V3 G ^. Let Ti := n S2, T2 := Si n S^, Ta := 5^ n := 
Si n S2, and Eij := {e = {x,y} G E \ x G Ti,y G Tj}. Using a sequence of estimates that 
is completely identical to (fT3])-(fT6l). we get 2 > J2eeS{Ti) + Yle€5(T3) ^e, again implying that 
™™{Z^ee<5(Ti) ^e, X^eg5(r.^) Xg} < 1, Contradicting the assumption that x satisfies all blossom in- 
equalities. □ 
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Theorem 15 For any even set of vertices in the plane, there is a fractional perfect matching x of 
minimum stabbing number that has an edge of weight at least 1/5. For any set of vertices in the 
plane, there is a fractional spanning tree x of minimum stabbing number that has an edge of weight 
more than 1/3. 

Proof. For both problems, consider a fractional vertex with a planar support graph. To see the claim 
for matchings, note that there must be a vertex with degree at most five; as the total weight for each 
vertex is 1, the claim follows. To see the claim for spanning trees, note that the total edge weight is 
n — 1, and the number of edges is at most 3n — 6, implying that the average weight is larger than 1/3. 
□ 

Theorem [TSlprovides the basic ingredient for an iterated rounding algorithm: at each iteration, fix 
the weight of an edge of maximum fractional weight to one, and re-solve the linear program. In each 
iteration, the number of edges with fractional weight is reduced, so we get an overall polynomial-time 
algorithm for finding an integral solution. 

Unfortunately, Jain's original proof only guarantees a constant-factor approximation for objective 
functions that arise as the (weighted) sum of the edge variables. However, the situation is different 
for our objective function which is a maximum over certain sums of edge variables, so an additional 
argument is needed for establishing a constant-factor guarantee. We are hopeful that this argument 
can be completed some time in the future ifTSl . As we demonstrate in our experimental paper |[T2l . the 
practical performance seems to be even better than the theoretically possible guarantees of 5 and 3. 

6 Notes and Conclusion 

We have presented the first algorithmic paper on stabbing numbers, resolving the long-standing open 
question of complexity, and providing an approach that appears to be useful in theory and in practice. 
There are a number of interesting open questions. 

Our proofs rely on a strong degeneracy of the vertex set, and it would be interesting to see a proof 
for vertices in general position. 

We were not able to extend our AAP-hardness proof to the case of finding a triangulation of 
minimum general crossing number; another interesting issue is how to convert the AAP-hardness proof 
for triangulations of minimum crossing number into an AAP -hardness proof for minimum stabbing 
number. 

Probably the most intriguing open question spawned by our work is whether the iterated rounding 
scheme suggested by the existence of a heavy edge in an optimal fractional solution to our linear 
programs (Lemmas [13] and [141 ) does indeed lead to a constant-factor approximation algorithm. Also, 
the use of the ellipsoid method (at least as a theoretical argument) is not "combinatorial," which always 
has to be considered a drawback. 

Another interesting question is to decide the existence of structures of small constant stabbing 
number. As the hardness proof for deciding the existence of a matching of stabbing number 5 illus- 
trates, this is still not an easy task. From some solvable special cases, we only note one; for a proof, 
see EH. 

Theorem 16 St-Tre2{P)=2 and St-Mat{P)=2 can be decided in polynomial time. 
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One may also ask for minimizing the average instead of the maximum stabbing number, and refer 
to the average over the whole continuum of lines intersecting a set of line segments, instead of just a 
combinatorial set of representatives. This, however, amounts to solving problems of minimum length, 
with all implications to hardness and approximation; again, see |[T2]| for details. 

Theorem 17 A set of line segments has minimum average (axis-parallel, resp.) stabbing number with 
respect to uniform distribution of lines, if and only if the overall Euclidean (Manhattan, resp.) length 
of all line segments is minimum. 

We remark that a linear program for minimizing the average stabbing number can be written with 
a sum in the objective function (instead of a maximum as we had to model it), allowing us to directly 
apply our iterated rounding technique and obtaining the desired approximation factors of 3 and 5, 
respectively; see iT2l . 
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